Constraining Interactions in Cosmology's Dark Sector 
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We consider the cosmological constraints on theories in which there exists a nontrivial couphng 
between the dark matter sector and the sector responsible for the acceleration of the universe, in 
light of the most recent supernovae, large scale structure and cosmic microwave background data. 
For a variety of models, we show that the strength of the coupling of dark matter to a quintessence 
field is constrained to be less than 7% of the coupling to gravity. We also show that long range 
interactions between fermionic dark matter particles mediated by a light scalar with a Yukawa 
coupling are constrained to be less than 5% of the strength of gravity at a distance scale of 10 Mpc. 
We show that all of the models we consider are quantum mechanically weakly coupled, and argue 
that some other models in the literature are ruled out by quantum mechanical strong coupling. 

PACS numbers: 



I. INTRODUCTION 

Multiple, complementary cosmological observations all 
suggest that the Universe has recently embarked upon 
an epoch of accelerative expansion. These observations 
include the cosmic microwave background (CMB), for 
example Refs. [H, 0, i, i, i, large scale structure 
surveys, for example Refs. [3,H, 0], including baryon 
acoustic oscillations and Type la supernovae 

[isl . [T^ . [TBI , [l^ . There now appears to be irrefutable 
evidence that the Universe's expansion deviates from 
that predicted by Einstein's General Relativity and a 
Universe solely populated by baryonic matter and radia- 
tion. 

Two new components, dark matter, that does not in- 
teract with light but does cluster under the force of grav- 
ity, and dark energy, that drives cosmic acceleration, have 
been invoked to resolve the disparities. In the minimal 
picture, dark matter does not feel any significant inter- 
actions, even with itself, apart from through gravity, and 
dark energy is a cosmological constant, not evolving and 
having no spatial fluctuations. Although this picture is 
wholly consistent with observations, the theoretical ori- 
gin of both of these dark additions still remains a mystery, 
and the simple interpretation above has its own issues, 
such as the coincidence and fine tuning cosmological con- 
stant problems. 

Recognizing that the physics of the dark sector is ef- 
fectively unknown at present, and in light of the pos- 
sible complexity of the dark sector arising out of high 
energy theory, theoretical models beyond the minimal 
picture have been considered. This includes a plethora 
of fundamental dark matter particle candidates, see 
for example (l7j for a review, that might well be ex- 



pected to have interactions beyond purely gravitational 
ones [ll,[li,[23,[2ll[23,[2l,[2l[2i,[li,[13 • Such interac- 
tions can have astrophysical consequences, for example, 
the prospect of dark matter interactions, such as self- 
annihilation, that could give rise to the 511 keV emission 
fH; the 'WMAP haze' m 



streams in galactic systems [32 



31| 



implications for tidal 
, as well as modifica- 
tions to dark matter halo profile [3j, [H, HI] , dark matter 
halo mass function [37j or altered dark matter motion in 
cluster collisions, such as the Bullet Cluster [sst . 

One possibility that could mitigate the cosmological 
constant problems is that non-minimal interactions ex- 
tend more broadly between dark sector particles, so that 
the properties of dark energy and dark matter arc cou- 
pled in some way. Such a direct coupling can be employed 
to address the coincidence problem, by relating the on- 
set of cosmic acceleration with the properties of a matter 
dominated universe [19, ,23, 39, 40, 41, 42, 43, 44, 45, 4§]. 
They can, however, also give rise to dynamical instabili- 
ties in the growth of structure [13, [H, Ei, [13, [HI, [Ell . 

The paper proceeds as follows: in section [ll] we de- 
scribe two examples of dark sector interactions, coupled 
dark matter-dark energy models in III Al and the Yukawa 
dark matter interaction in III Bl that can have astrophysi- 
cally observable consequences. In III CI we summarize the 
theoretical and observational constraints on dark sector 
interactions, a subset of which we focus in on detail in the 
paper. We present the constraints from the latest cosmo- 
logical observations on coupled dark matter-dark energy 
models in section IIIII and the Yukawa dark matter in- 
teraction in IIVI In section |V] we discuss the restrictions 
placed on models in the strong coupling regime. Finally, 
we pull together our findings and discuss their implica- 
tions in section IVll 



II. INTERACTING DARK MATTER 
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In this paper we consider scenarios in which a purely 
dark sector interaction exists, resulting from a non- 
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minimal coupling of dark matter to a scalar field. Such 
couplings give rise to additional forces on dark matter 
particles in addition to gravity. In this section we de- 
scribe two examples of models that exhibit this behavior. 
In the following sections we will discuss the observational 
constraints on these models. 



A. Coupling dark matter to dark energy 



Consider the general action 




1/(0) 



(2.1) 



where g^i, is the metric, Mp = {87rG)~^^'^ is the reduced 
Planck mass, and we use natural units with h = c = 1. 
Here </> is a scalar field which acts as dark energy, "ifj are 
the matter fields in the jth sector described by the action 
Sj, and aj{4>) describes the coupling of the scalar field to 
the jth sector. This general action (|2.ip describes a wide 
range of models, inclu ding the Einstein frame version of 
f{R) modified gravity [53l. Is^. Issl. IsgI. ISTj . A special case 
is when the couplings are identical in all the different 
sectors, OLj{(j)) = a{(f>) for all j, in which case the theory 
satisfies the weak equivalence principle. 

Although violations of the equivalence principle are 
strongly observationally constrained for normal matter, 
the constraints on dark matter are much weaker, as 
emphasized by Damour, Gibbons and Gundlach [2]| . 
Therefore it is interesting to consider models with two 
sectors, dark matter with coupling function ac((/>), and 
normal (baryonic) matter with coupling function ah{4>). 
Such models will automatically satisfy observational con- 
straints on the weak equivalence principle that involve 
only baryonic matter. They must also satisfy the addi- 
tional constraint from Solar System observations that 



(2.2) 



where 0o is the present day cosmological background 
value of ^. Below we will specialize to models with 
ttf, = 0, in which the scalar field is coupled only to the 
dark matter, which automatically satisfy the solar system 
constraint (|2.2p . 

We note that theories of the form (|2.1[) , in which differ- 
ent sectors couple in different ways to the scalar field 
arise very naturally from higher dimensional models with 
branes. An example is provided by the Randall Sundrum 
I (RSI) model [63| . with two parallel branes in a five di- 
mensional anti-deSitter space, one with positive tension 



^ This assumes that the solar perturbation to (f) is in the Unear 
reg ime, which is not true, for example, in chameleon models 
[5l|5i,|63,[63,[6l. 



and one with negative tension. The low energy four di- 
mensional description of this model is of the form (|2.ip 
with no potential [6^ . [65l | , with two sectors correspond- 
ing to matter on the two different branes, which we will 
denote -I- and — . In this case the scalar field is a radion 
field that encodes the distance between the two branes 
in the fifth dimension. The two coupling functions are 

a+(0) lncosh(0/V6Mp), (2.3a) 
a-{4>) = lnsinh(0/\/6Afp) . (2.3b) 

The conventional interpretation of this RSI model is 
that visible matter lives on the negative tension branc, 
and that the positive tension ("Planck") brane contains 
a hidden sector. This interpretation requires that the ra- 
dion be stabilized, otherwise the Solar System constraint 
(|2.2p is violated for all values of the present day cosmo- 
logical value 00 of the scalar field. An alternative in- 
terpretation (which unlike the conventional one does not 
solve the hierarchy problem) is that visible matter is on 
the positive tension brane and that dark matter is on the 
negative tension brane, i.e., we make the identifications 
+ = b and — = c, and the radion is not stabilized. In this 
model, normal matter is minimally coupled in the limit 
of small 00 (corresponding to distant branes), so that the 
constraint (|2.2p can be satisfied in that regime. 

In the remainder of this paper we assume zero 
baryonic-scalar coupling, af, = 0, and we will denote the 
dark matter coupling function ac(0) simply as a(0). The 
model will then be specified completely by a choice of 
coupling function a(0) and potential V{(j)). 



B. Yukawa interaction between dark matter 
particles 

Rather than coupling dark matter to dark energy, we 
can also modify the coupling of dark matter particles 
with themselves. One class of models of this type involve 
an interaction between fermionic dark matter, ij), and an 
ultra-light pscudo scalar boson, 0, that interacts with the 
dark matter through a Yukawa coupling with strength 
described by the Lagrangian [33 |. 



+,90'0-0. 
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(2.4) 



For g 7^ 0, on scales smaller than ~ "i^^, the Yukawa 
interaction acts like a long-range 'fifth' force in addition 
to gravity. The effective potential felt between two dark 
matter particles is 



V{r) 



Gr 



1 + QfYuk exp 



with 



avuk 



(2.5) 



(2.6) 
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In our investigations of this model in Sec. IIVI we will 
neglect the cosmological effects of the scalar field (j>, and 
assume that dark energy is a cosmological constant^. The 
cosmological implications of Yukawa-like interactions of 
dark matter particles have previously been considered 
across a range of astrophysical scales, including dark mat- 
ter halos [12 m, [G^I , tidal tails [H, [sl] , cluster dynamics 
[sst , and large scale structure surveys [6l| . 

C. Theoretical and observational constraints 

Models such as the ones described above face a range 
of theoretical and observational constraints arising from 
both particle physics and gravity. We will focus on a 
subclass of these in this paper, but it is worth mentioning 
the general web of desiderata and constraints. These 
include: 

• The existence of an ultraviolet (UV) completion. 
Ideally one would like to find an embedding of the 
theory (|2.ip in string theory. Such embeddings 
have been recently found for inflationary models, 
see, for example, the review [g^. However it is 
difficult to find UV completions for quintessence 
models; see, for example, the supergravity no-go 
theorem in Ref. [tO]. 

• Fine tuning and the taming of loop corrections. 
Typically one would like a dark energy model to 
provide the unnaturally small value of the vacuum 
energy today. Having chosen such a small param- 
eter value in one's Lagrangian. it is often neces- 
sary to fine tune the model to prevent renormal- 
ization of parameters through couplings to other 
fields. This is sometimes avoided in dark energy 
models by making the dark energy field a pseudo- 
Nambu-Goldstone boson, such as in the Yukawa 
scenario discussed in III Bl This is not necessarily 
the case for the action (|2.ip . In this paper we shall 
just assume that such tunings exist in (|2.ip . since 
avoiding them is not our focus. 

• The strong coupling problem. If we treat the La- 
grangian (|2.ip as an effective field theory (as we 
should), valid up to some energy scale A, then there 
will exist irrelevant operators suppressed by pow- 
ers of the cutoff. In certain regimes, these opera- 
tors may become important, meaning that we are 



^ We note that the action 112.41 1 is actually a specific case of our 
general action 1 12.11 1. with V{(l>) = rrv^cjp /2 and a{4>) = ln[l — 
gij>/m^\/'i, and specialized to the regime where the fermions are 
non-relativistic so that we can neglect the modifications to the 
fermion kinetic term in the action. However, our interpretation 
of this model is different from our interpretation of the models 
1)3. 9|l and l|3.16|l . since the (fi field is not the dark energy and we 
are neglecting its cosmological evolution. 



no longer able to trust the effective theory. This 
will not arise in the theories we discuss here in a 
cosmological context. This strong coupling issue is 
discussed below in Sec. |Vl 

• Disagreement with the required background cosmol- 
ogy. Obviously, a successful model must be able 
to reproduce the correct expansion history of the 
universe, preferably without excessive fine tuning 
of initial conditions. This can be a real problem 
for some models, for example some f{R) modified 
gravity models [71[ . In sections IIIII and IIVI we in- 
vestigate cosmological evolution in coupled models. 

• Problems with linear perturbations around the FRW 
solution. Here the possibilities include disagree- 
ments with solar system tests of gravity [56|, or 
incorrect predictions for the linear power spectrum 
of matter perturbations. In addition instabilities 
causing catastrophic collapse of over-densities can 
be pre sent in some regimes for coupled theories 

113, 11, Si, Sim, [13]. 

• Problems in the nonlinear regime There is the also 
possibility of interesting phenomena in the nonlin- 
ear regime. Some may be positive; for example the 
Chameleon effect jssl . [59l | can ameliorate problems 
with Solar System tests (t^, [ll] ■ Some other phe- 
nomena can be problematic, for example in some 
models the spatially averaged metric is not a so- 
lution of the field equations that one obtains by 
assuming homogeneity and isotropy (i.e. the "mi- 
croscopic" and "macroscopic" field equations differ) 

Bill. 

In this paper we will focus on the constraints obtained 
from the background cosmological evolution, linearized 
cosmological perturbations, and the strong coupling con- 
straint. 



III. COSMOLOGICAL CONSTRAINTS ON 
COUPLINGS BETWEEN DARK MATTER AND 
DARK ENERGY 

In this section we consider the class of models 
specialized to two sectors, the visible sector with zero 
coupling function, and the dark matter sector with cou- 
pling a(0) [n, [ii, SO, sn, m m Hi, H, s [5i|. The 

resulting equations of motion are 

-"^'('/')3a6 + e"(*Vc"a^ih, (3.1a) 
VaV>-r((/)) = a'(0)e"('^Vc, (3.1b) 

^aiPcU'') - 0, (3.1c) 

u^X/bu" = -a'(</.)(g"^ + w"M'')V6(/). (3. Id) 

Here Gab is the Einstein tensor. Tab is the stress-energy 
tensor of visible matter and m° is the four velocity of 
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the dark matter. The quantity pc is proportional to the 
number density of dark matter particles with respect to 
the metric gab, it scales cx like uncoupled dark matter 
in the background cosmological solution. The observed 
energy density of dark matter is e"pc- 



A. Evolution of background cosmology 

Writing the flat FRW metric as 

ds^ = a'^{T){-dT^ +dx^) , 

with scale factor a(r) and conformal time r, the Fricd- 
mann equation is 

SM^H' = i02^aV(0)+a2e"(*Vc+aVh+aVr , (3.2) 

where dots represent derivatives with respect to r and 
Ti. = a I a. Here and pr are the densities of baryons 
and radiation. The remaining equations for the system 
are 

+ 27i:(/) + aV(0) = -a2a'(<?!))e"('*Vc, (3.3a) 
Pc + iUpc = 0, (3.3b) 
Pb + "iTipb = 0, (3.3c) 
Pr + mpr = 0, (3.3d) 



dropped, yields 

dx 

dN 

dy 

dN 
dX 
dN 



-\/6 ,93, 9 

-3x + —Ay + -x(l + x- 

+C{1 - .t2 - y^) , 

V6x 3 , 9 2\ 
— Y'^xy+ -y{l + x -y ) 



(3.6a) 
(3.6b) 

(3.6c) 



In these equations, F and C are understood to be the 
functions of A obtained by eliminating (f> in Eqs. (j3.4p and 
p.Sp . The fixed points of this system are the solutions 
of the equations dx/dN = dy/dN = dX/dN = 0. 

After Eqs. (|3.6p have been solved to obtain the func- 
tions x{N), y{N) and A(7V), the Hubble parameter HiN) 
can be found from Eqs. (|3.4p . and the dark matter den- 
sity Pc can be obtained from the Friedmann equation, 



3A/27^2 



1 . 



(3.7) 



Note that the effective total equation of state parameter 
Woff, defined by dlna/dlnr = 2/(1-1- 3wes), is simply 
given by 



where primes denote derivatives with respect to 



WcS 



(3.8) 



1. Dynamical attractors in general coupled models 

Scalar field quintessence models of dark energy have 
been shown to have expansion histories that exhibit scal- 
ing attractor solutions which reduce sensitivity to initial 
conditions for the scalar field [H, S S H [Z^. The 
same has been found to be true of coupled quintessence 
scenarios [H, 0, EH, EH, f{R) gravity [ni and scalar- 
tensor gravity [80[. 

We specialize to the matter dominated era and neglect 
the baryons and radiation. To describe the attractor be- 
havior in coupled models described by Eqs. p.2p - (|3.3p . 
we use the dimensionless variables defined by Copeland 
et. al. [si, El: 



y 



aVV 



X: 



' V ' 



VV" 
(3.4) 



along with the dimensionless coupling variable, 

C{4>) ^ (3.5) 



/3 



with [3 = ■\/2/3. Rewriting the evolution equations 
- p.3ap in terms of these variables and in terms of the de- 
pendent variable = ln(a), with baryons and radiation 



from Eqs. (EH) - (IHTTl) . 

We consider the dynamical behavior for two specific 
models, with an exponential and power law potential, 
in the presence of an exponential coupling between the 
scalar field and cold dark matter. 



2. Model 1: An exponential potential 

We will consider a model with an exponential potential 
and linear coupling given by 



C(3<j> 
Mr, ■ 



(3.9a) 
(3.9b) 



Here A and C are dimensionless constants of order unity, 
and Vo is a constant of order M^Hq. For this model the 
functions r(0), A(0) and C{(j)) defined by Eqs. and 
p.Sp are constants: 



r(0) = 1, A(0) = A, C(0) = c. 



(3.10) 
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There exist three fixed points (these are a, 6m, Cm from 
Amcndola's analysis of this specific model [4l|): 



{x,y) 



2C ■ 



V6' \ 6 



(3.11a) 
(3.11b) 



1/2 



1 



X-I3C ' \2j A-/3C 



1/2 



X (1 + (3*C^ - /?3CA) 



1/2 



(3.11c) 



It is important to note that, depending on the values of 
the parameters of the model, some of these fixed points 
may not exist, i.e., they may be complex rather than real. 
In addition, when they do exist, they may or may not be 
stable attractors during the matter and dark energy era. 

The first of these (|3.11ap is an attractor approached 
as the matter dominated era is entered. The potential is 
subdominant and the scalar field kinates, leading to an 
effective equation of state parameter 



Weft 



4C^ 
~9~' 



(3.12) 



3. Model 2: A power law potential 
We also consider the inverse power law potential model 



F(0) 
a(</>) 



Foexp(^— 



where n is a constant for which 



r(0) = ^, Xi4>) = -n('-^ 



(3.16a) 
(3.16b) 

C{4>) = C. 



(3.17) 

There are two stable attractors which arise in the matter 
and accelerated eras, respectively. 



2C 



{x,y) = [0,1], 



(3.18) 
(3.19) 



for which, in both cases, A ^ 0. Eq. (|3.18p gives a matter 
dominated era attractor equivalent to p.l2p , while p.l9p 
is an accelerative attractor with Wcff = — 1, independent 
of C and n. 



This evolution is often described as a '0CDM' era, and 
differs from the usual CDM dominated era with WeS = 0. 
Its existence and properties have led to significant issues 
when fitting some f{R) theories, for which C = 1/2, to 
observations [7l|. 

The second fixed point, (|3.11bp . is a stable attractor 
for 

A(A - /3C) < 3 (3.13) 
with effective equation of state parameter 

A^ 

w;eff = -l + y. (3.14) 

This attractor gives rise to acceleration if A^ < 2. This 
fixed point arises entirely from the nature of the scalar 
potential, and is independent of the coupling C; in par- 
ticular it arises in the minimally coupled case C = 0. 

The final fixed point, (|3.11cp . with 

exists if the second fixed point is unstable. We will find, 
however, that condition (j3.13p is satisfied in the viable 
models we analyze below, so that this final fixed point 
does not arise. 



4- Numerical evolution of attractors 

In Figure [T] we show the background expansion history 
for examples of the exponential and power law potentials 
and the coupling discussed here. 

Typically in these models, the radiation era evolution 
is the same as in ACDM, with scalar field attractors with 
fi^ = or Wcj, = 1/3. In certain cases, e.g. exponential 
models with A > 2, however, the radiation era can be re- 
placed by a kinetic scalar field dominated era for models 
with Hq consistent with HST. However these models do 
not confront data well. 

A difference between ACDM and the coupled scenarios 
can arise in the matter dominated era as described above. 
In this regime the attractor evolution alters the matter 
dominated expansion history via equation (|3.12p . The 
angular diameter distance of the CMB and the growth 
functions for large scale matter perturbations (A; < kgg) 
entering the horizon after matter radiation equality, rel- 
ative to the smaller scale (k > k^q) perturbations, is al- 
tered in comparison to ACDM. 

At late times, the coupled models tend towards accel- 
erative attractors which are independent of the coupling 
C, given by (|3.11bp for the exponential and p.l9p for the 
power law potentials, respectively. Note, however, that 
the evolution will not necessarily have reached the attrac- 
tor today, and the coupling can therefore play a role in 
determining ifoff by altering the time at which the shift 
from the (pGDM to the accelerative attractor occurs. 
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V(0)~exp(-0/Mp) V(0)~V0 




FIG. 1: Examples of evolution of the effective equation of state, WcS, in coupled scalar field dark matter models with an 
exponential potential V{(j)) oc exp(— ^/Mp) (left panel) and a power law potential V{<l)) oc 1/0 (right panel). Cosmological 
parameters are fixed to Ho = 70, ilc = 0.25, fit = 0.05, and C = 0.1 (black) and C = 0.5 (red). Both models follow the 
coupling dependent attractor in the matter dominated era and asymptote to coupling independent attractors at late times. 
The timing of the transition between these two attractors is sensitive to both the potential and coupling parameters. For 
the exponential potential the dynamical attractor leads to a negligible dependence on initial conditions, shown here through 
comparing evolution with two different initial values of 0i = <j}{a = 10~*), = lAfp (full) and 10~^°Mp (dashed). For the 
power law potential, however, a sensitivity to initial conditions can exist in the transition era. This is accounted for in the 
analysis by marginalizing over initial conditions. 



As shown in Fig. [1] for the exponential potential the 
attractor behavior quickly takes over, and the initial con- 
ditions have no effect on the dynamical evolution. In the 
case of the power law potential, however, we find there 
can still remain some sensitivity to the initial value of the 
scalar field during the transition between matter domi- 
nated and accclcrative attractors. As discussed in section 
IIII C[ we account for this in the analysis by marginalizing 
over the initial value of (h. 



B. Evolution of linearized cosmological 
perturbations 

As well as background evolution, we are also interested 
in the predicted evolution of density perturbations. We 
write the inhomogeneous density and scalar field as 

Pcix,T) = Pc{t){1+6c{x,t)), (3.20a) 
(t){x,T) = (j){T) +ip{x,T). (3.20b) 

We use the notation of Ref. [8l| to describe the perturbed 
metric in synchronous gauge in terms of two functions 
77(t) and /i(t). The four independent components of the 
Einstein equation are then 

2k^r]-nh = -a^e°'pc{S + a'ip)~a'^V'ip 

-(j)(p, (3.21a) 
2k^7) = a^e^pcOc + k^'ixp, (3.21b) 
h + 2nh-2k^r] = -3^ip + 3a^V'(p, (3.21c) 



and 



6ii + h + 2H{h + 6t]) ~ 2k'^r] = 0. (3.22) 



Here k is the comoving wavevcctor and 6c ~ ikjV^^ is the 
gradient of the CDM peculiar velocity, Vc- Also we have 
specialized to units with Mp — 1. We include just the 
effects of CDM and the scalar field, and neglect baryons 
and radiation, for simplicity. The perturbed fluid equa- 
tions are 



Sc + h + 9c = 0, 



dT 



(3.23) 
(3.24) 



'f + 2H^ + [k^ +a^V" + a'^e°'pc{a" + {a'f)]if 
= a^a'e°'ScPc- (3.25) 

There exists an extra gauge degree of freedom that 
preserves synchronous gauge, given by the the coordinate 
transformations 



(3.26) 



dT 



kcoTZlikje'"'-''] / — , (3.27) 



where Cq is a constant and k^. = kj/k. Under this trans- 
formation the metric and matter perturbations transform 
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FIG. 2: Joint 68% (dark shaded) and 95% (light shaded) constraints in the exponential potential model for the the fractional 
matter density, fim and the effective scalar equation of state, w^{leit panel), and the coupling, C (right panel). The comple- 
mentary constraints arising separately from the WMAP CMB spectra, SDSS matter power spectrum plus SDSS and 2dFGRS 
baryon acoustic oscillation data sets, and the 'union' Type la supernovae data set, and a HST prior on Ho are shown. 
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We can define two new variables 



5c + 3H- 



(3.28) 
(3.29) 

(3.30) 

(3.31) 
(3.32) 

(3.33a) 
(3.33b) 



that are invariant under the residual gauge transforma- 
tions. 

In the minimally coupled case «(</>) = 0, one typically 
fixes the residual gauge freedom by choosing the CDM 
rest frame in which 0c = and 5c = [Sil- This is 
consistent since, once fixed to zero, 9c remains zero at 
all times, by Eq. (|3.24p . In the presence of an evolv- 
ing non-minimal coupling, however, if the CDM velocity 
divergence is initially zero it will evolve to be non-zero. 
Therefore it is not possible to fix the residual gauge free- 
dom in this way. In our computations below we will 
evolve the perturbation equations in an arbitrary syn- 
chronous gauge, and then use Eqs. (|3.33|) to pick out 
gauge invariant combinations of the perturbation vari- 
ables. 



C. Comparison with data 

We have modified the CAMB code [111 to evolve 
background equations and first order density perturba- 
tions for a flat universe containing baryons, CDM. radi- 
ation, massless neutrinos and a scalar field coupled non- 
minimally to CDM and use CosmoMC [13] to perform a 
Monte Carlo Markov chain of the model parameter space 
in comparison to current cosmological data. We explore 
the exponential model p.9p and the power law model 
p.l6p allowing the exponent A and index n to vary in 
each case. We set the initial conditions for the scalar 
field well into the radiation era, a = 10~^, allowing the 
initial value of the scalar field to vary, then evolve for- 
ward to the times at which CAMB typically begins its 
integration. 

We constrain the models using a combination of cos- 
mological datasets, including the measurements of the 
CMB temperature and polarization power spectrum from 
the WMAP 5-year data release [1, [s^, the 'union' set of 
supernovae compiled by the Supernovae Legacy Survey 
(SNLS) [11] , and we impose a Gaussian prior on the Hub- 
ble constant today, Hq = 72± 8, using the Hubble Space 
Telescope (HST) measurements [13]. We use the mat- 
ter power spectrum of Luminous Red Galaxies (LRG) as 
measured by the Sloan Digital Sky Survey (SDSS) survey 
[1,13, for which we include the shift parameter, a^d, to 
adjust the matter power spectrum as discussed in 



a-sci 



dv 



dv{z = 0.35)('"°'^'='' 

dv{z = Q.i^)U^duc^al) 

{l + zfdA{zfcz^^'^ 



H{z) 



(3.34) 



(3.35) 



where dA[z) is the physical angular diameter distance at 
a redshift z and the fiducial model is a standard ACDM 
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FIG. 3: Joint 68% (dark shaded) and 95% (light shaded) constraints using combined WMAP CMB, SDSS matter power 
spectrum, SDSS and 2dFGRS baryon acoustic oscillation, 'union' type la supernovae datasets and HST Ho prior for the power 
law potential (blue) and exponential potential (red) for the the fractional matter density, fim and the effective scalar equation 
of state, ui0, (left panel), and the coupling, C (right panel). 



model with il™ — 0.25, J^a — 0.75 and with the same 
Hubble constant as the theory model. We also use con- 
straints on the expansion history from the Baryon Acous- 
tic Oscillation data of the 2dFRGS and SDSS surveys 
based on measurements of the ratio of the sound 
horizon at last scattering, rs{z^), to the distance mea- 
sure dviz) at z = 0.2 and z = 0.35. Since the dynamical 
attractor solutions, in the presence of a non-minimal cou- 
pling, can alter the background evolution in the matter 
dominated era, one finds that the redshift of last scat- 
tering, z*, can no longer be accurately estimated using 
the fitting formula of Hu and Sugiyama [s^. Instead 
we calculate the redshift of maximum visibility and use 
this as the appropriate measure for the redshift of last 
scattering. 

In Figure [2] we show the complementary 2D marginal- 
ized constraints for the exponential potential model in 
light of the various cosmological datasets. The CMB 
data (along with the HST prior on Hq) provide the best 
individual constraint on the coupling strength with ID 
marginalized constraints |C| < 0.13 at the 95% confi- 
dence level (c.l.). The Type la supernovae alone provide 
only weak constraints on both the coupling and on the 
total matter density in a non-minimally coupled model. 
This is because the coupling allows a late time, cosmo- 
logically consistent expansion with WcS ~ ^<fiW^ ^ —0.7 
to be generated by a strongly phantom-like model, with 
W0 <C — 1 and fi^ sa 0, where 



W4. 



02 _ V{C^) 



1^2 



102 



102 



(3.36) 



(3.37) 



V{(t)) + e"<^pc + Pb + P-i 

These models arc not consistent with CMB and LSS ob- 
servations, however. 



In Figure [3] we show 2D marginalized constraints for 
the exponential and power law potential models from all 
the cosmological datasets combined. The ID marginal- 
ized constraint on the coupling in the exponential poten- 
tial case is \C\ < 0.037(0.067) at the 68% (95%) con- 
fidence level. This represents a significant tightening of 
constraints over previous analyses, for example [sH ] found 
|C| < 0.1 at the 98.6% level using CMB data from the 
Boomerang satellite. The potential exponent. A, is con- 
strained to be |A| < 0.95 at the 95% c.l.. 

In the power law potential case the ID marginal- 
ized constraint on the coupling is comparable, with 
-0.026(-0.055) < C < 0.034(0.066) at the 68% (95%) 
C.L. Again, the constraints on the coupling strength 
have improved with the increased precision and comple- 
mentary variety of the cosmological data, e.g. a pre- 
vious analysis with first year WMAP data alone found 
C < 0.085(0.159) at the 68% (95%) c.l. Within the 
range investigated, —6 < n < 6, the power law exponent, 
n, is not constrained by the data. 

In both the exponential and power law potential cases 
the constraints are wholly consistent with a minimally 
coupled ACDM model (A = or 7i = 0, and C = 0) at 
the la level. 



IV. COSMOLOGICAL CONSTRAINTS ON A 
YUKAWA-TYPE DARK MATTER 
INTERACTION 

The astrophysical implications of Yukawa-like interac- 
tions have been considered across a range of scales: in 
the context of dark matter halos [H, [6^, [fi^ ; tidal tails 
[H, nil ; cluster dynamics [111 ; and large scale structure 
surveys [68j . In our analysis we consider large scale cos- 
mological constraints on a Yukawa coupling described in 
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FIG. 4: CMB temperature power spectrum comparing a fiducial minimally coupled ACDM model (grey) with four models 
with the Yukawa interaction with oyuIc = 10, rB=2Mpc (red), ayuk = 5, rs = 2 Mpc (green) avuk = 10, r-i,=5Mpc (magenta), 
CfYuk = 5,rs=5Mpc (blue). Data from WMAP 5 year (blue points) and ACBAR (black points) experiments are also shown. 
The inset plot shows a blow up of the small scale anisotropics measured by ACBAR. 



section I IIBI We modify the publicly available CAMB 
code [Sst to include this modified force between dark mat- 
ter particles. This alters the growth of matter perturba- 
tions. For example, the dark matter density fluctuations 
evolve according to 



Gcik) 



G 



PcSc + PbSb + 2pjS^ 



0. (4.1) 



Here Gc{k) is the effective gravitational constant govern- 
ing the interaction between dark matter particles, given 
from Eq. ([231) by 



Gcik) = G 



1 



avuk 



1 + (krs 



(4.2) 



We use the CosmoMC code [S^ to obtain cosmological 
constraints on the ratio Gc/G from the 5 year WMAP 
CMB temperature and polarization data p, |83|, small 
scale CMB temperature data from ACBAR [j, Q and 
the SDSS LRG matter power spectrum [1]. We include 
CMB lensing, and marginalize over the amplitude of the 
secondary Sunayev-Zel'dovich anisotropics. 

In Figure S] we show the effect of the Yukawa cou- 
pling on the CMB temperature anisotropics. With the 
addition of small scale anisotropy measurements from 
ACBAR, constraints on the interaction are able to be 
made. 



In Figure [H] we show the constraints on Gc/G at two 
scales, 1 Mpc and 10 Mpc with Gc/G{lMpc) < 2.7 and 
Gc/G{10Mpc) < 1.05 at the 68% confidence hmit. The 
improvement in the fit to the data obtained by introduc- 
ing the Yukawa interaction is not statistically significant 
however, the best fit effective = — 21n£ = 1354.0 in 
comparison to 1354.1 for a ACDM model. 

Yukawa interactions on the levels allowed by large 
scale constraints could well have interesting implications 
for gravitational dynamics on cluster, galactic and sub- 
galactic scales [12, m, 113, HEl- Frieman and Gradwohl 
[34| argue that the intracluster gas distribution could 
constrain — 0.5 < avuk ^1-3 for of a few hundred kpc, 
which would translate to -0.5 < Gc/G'(lMpc) < 2.2, 
comparable with our constraints from large scale data. 
Kesdcn and Kamionkowski (3^.[33j demonstrate that cou- 
plings of strength Gc/G > 1.04 on < 100 kpc scales could 
well have observable implications for baryonic and dark 
matter distributions in tidal disruptions of dwarf galax- 
ies, although a comparison with data is yet to be per- 
formed. We leave a detailed analysis of the joint con- 
straints on Yukawa interactions from combined astro- 
physical and cosmological scales to future work. 

We note that the observational constraints on the 
Yukawa coupling ayuk also yield constraints on the more 
general class of models (|2.1[) discussed in Sec. IIIIl pa- 
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FIG. 5: ID likelihood constraints on Gc/G at IMpc (left panel) and lOMpc (right panel) for the Yukawa dark matter interaction, 
in light of WMAP 5 year and ACBAR CMB anisotropy, and SDSS LRG matter power spectrum observations. 



rameterized by a baryonic coupling function ab{4>) and 
a dark matter coupling function ac{4>). In these models 
the effective Newton's constant G^- for coupling between 
sector i and sector j is given by Gij = G{1 + ji^j) with 
7i — ^/2Mpa^{(f>o) [Hlj. Now dark matter is observed only 
through its gravitational interactions. Therefore the ob- 
servations cannot distinguish between a situation with 
baryonic and dark matter densities pb, Pc and Newton's 
constants Gcc, Gcb and Gbb, and a situation with den- 
sities Pb, Pc and coupling constants e~'^'^Gcc, e^^Gcb 
and Gbb, where ly is an arbitrary constant ^. If we define 
the parameter 



GrrG, 



a 



bb 



^cb 



(l + 7g)(l + 7g) 
(1 + 7c7fc)' 



-1, 



(4.3) 



then we see that a is invariant under the above symmetry, 
and also a reduces to ayuk for the models discussed in 
this section for which = 0, at short lengthscales r ^ 
Ts- It follows that the arguments of Ref. [13] give the 
constraint 



0.5 < a < 1.3 



(4.4) 



on the class of models (|2.ip . This constraint already sig- 
nificantly limits some models, for example together with 
the Solar System constraint (|2.2p it rules out the version 
of the RSI model we discussed in Sec. Ill Al above. 



QUANTUM MECHANICAL STRONG 
COUPLING CONSTRAINT 



General relativity is a weakly coupled effective quan- 
tum field theory at lengthscales large compared to the 
Planck length [sol. [90j. However, many modifications of 
general relativity do not share this property. It can hap- 
pen that at relatively low energies, loop corrections be- 
come large and one can no longer trust the classical the- 
ory. The theory becomes strongly coupled, like quantum 
chromodynamics at low energies. This occurs for theories 
of massive gravitons at energy scales above (m^Afp)^/^, 
where nig is the graviton mass and Mp is the Planck 
mass [9l|, and in the DGP model at lengthscales below 
~ 1000 km 



92l |. It is also generic for theories which 



modify gravity in the infrared without introducing new 
degrees of freedom . 

Many coupled cosmic acceleration models in the litera- 
ture are invalid because of this consideration. There is a 
straightforward procedure for computing when a model 
of the form (|2.ip is in the strong coupling regime: for 
a given classical solution, compute the action of fluctua- 
tions S(j) about that classical solution, and then Taylor ex- 
pand that action. The Taylor expansion of the potential 
gives nonrenormalizable terms of the form (50"+^/ (A„)", 
where A„ is some energy scale and n > 1 is an in- 
teger. Then the theory is strongly coupled at energy 
scales above the lowest of the scales A„, ie for E > Age 
where Age = min„A„. If the corresponding lengthscale 
r ^ 1/Asc is in the range probed by observations, then 
the predictions of the classical theory cannot be used to 
compare with observations. We will see that some models 
are ruled out by this consideration. 

For a generic scalar field theory which acts as a model 
for cosmic acceleration, the potential can be written as 



^ Note that Gcc is denoted Gc in the rest of this paper. 



y(0) ^ HlMlV{cp/Mp) 



(5.1) 
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where the function F is a dimensionless function of a di- 
mensionless argument. For a generic model we expect 
the derivatives of V to be of order unity, so the nth 
term in the Taylor expansion scales as Hq {(j) / Mp)^ . 
The corresponding strong coupling scale A is A ~ 
Mp^Mp/Hny^^^-"-^ which is larger than Mp for n > 4. 
Thus for generic quintessence models there is no strong 
coupling issue just from the potential. 

We next discuss the effects of matter coupling. Con- 
sider a generic theory of the form (|2.ip with coupling 
function a(0) and potential V{(f>), for which the poten- 
tial contains a nonrcnormalizable term 



(5.2) 



with 71 > 1. Wc first show that, for a localized source of 
mass ^ M and size ~ R, such a term has a significant ef- 
fect classically before it leads to strong coupling, as long 
as the mass M is sufficiently large. Thus, in the regime 
where one can treat perturbations from matter inhomo- 
geneities linearly, the theory is never strongly coupled. 

To see this, we denote by 6(f> the perturbation to the 
cosmological background solution (j){t) that is generated 
by the localized source. We take the ratio of the terms 
{S(t>)^+'"/A" and {VScf))^ in the action to get 



A" 



(5.3) 



We assume that in the absence of the term (|5.2p in the 
action, the scalar field can be treated as a massless field 
with dimensionless coupling strength C to matter, obey- 
ing an equation of the form 



OS(t> = -/3Cp/Mp. 



(5.4) 



This gives the following order of magnitude estimate 
(dropping factors of order unity, including /3) for the 
value of the scalar field perturbation near the surface at 

r^R, 



\C\M 
MpR ■ 



(5.5) 



The ratio (|5.3p will therefore be small for A ^ Ad, where 
the critical value of A for the perturbation to the poten- 
tial to be important classically is 



A, 



1_ f \C\M 
R \ M., 



l+2/r, 



(5.6) 



Next, the theory will be strongly coupled at energies 
E > A, and for a source of size ^ R the relevant quanta 
have energies E ^ 1/ R, so strong coupling will occur for 
A > Asc = l/R. Comparing this with the estimate ([57 
we see that as long as 



\C\M > Mj 



(5.7) 



the theory will never be strongly coupled in the regime 
where the perturbation to the potential can be neglected 



classically. The condition (|5.7p will be satisfied for all 
astrophysically relevant sources when C is of order unity. 

This analysis applies to a large class of scalar field the- 
ories, including the models discussed in this paper. How- 
ever it does not apply to theories where nonlinearities are 
important classically, such as chameleon field models and 
f{R) versions of these. The strong coupling constraint 
on these models must be checked on a case by case basis. 
For example, for the /(i?) theory of Faulkner et al. [t^ . 
an order of magnitude estimate gives the strong coupling 
scale to be A^c ~ M(M3Afp/p)i/("+i) where p is the 
matter density and M , n are parameters of the model. 
Demanding that Age not be so low as to invalidate the 
predictions of the model rules out a large portion of the 
parameter space in this case. 



VI. CONCLUSIONS 

The cosmological observations of the past few decades 
have provided firm evidence for significant physics be- 
yond the standard model of particle physics. It now 
seems clear that the successful formation of structure in 
the universe demands a new particulate component of 
the cosmic energy budget - dark matter - and that cos- 
mic acceleration may require some kind of dark energy, 
or a significant infrared modification of general relativity. 

While these phenomena have been revealed through 
their gravitational effects, their microphysical properties 
remain undetermined although, of course, there exist 
many complementary bounds on what those properties 
may eventually prove to be. A priori, there is no rea- 
son that dark matter and the physics responsible for cos- 
mic acceleration should find themselves in entirely dis- 
connected sectors of our underlying fundamental theory. 
Indeed; that small portion of the energy budget about 
which we know a great deal - visible matter - comprises 
a richly detailed spectrum with multiple interactions and 
a beautiful underlying symmetry structure - the standard 
model. It thus seems reasonable, in light of our current 
ignorance regarding the nature of cosmology's dark sec- 
tor, to explore possible interactions between dark matter, 
cosmic acceleration and visible matter. 

In this paper we have studied a broad class of cou- 
pled dark matter dark energy models, and have inves- 
tigated the constraints on such models from the most 
recent supernovae data, from precise measurements of 
the large scale structure of the universe, and from cos- 
mic microwave background experiments, including the 
recent WMAP 5-year data release. While it is notable 
how constraining each of these sources is individually, the 
combined constraints are surprisingly strict. Indeed, for 
the class of models studied here, we have demonstrated 
that the strength of the coupling of dark matter to a 
quintessence field is constrained to be less than 7% of 
the coupling to gravity at the 95% confidence level. 

Furthermore, we have applied our techniques to mod- 
els of possible astrophysical interest, in which long range 
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interactions between fermionic dark matter is mediated 
by a light scalar with a Yukawa coupling. We have 
shown that large scale cosmological measurements con- 
strain such interactions to be less than 5% of gravity at 
a distance scale of 10 Mpc. 

We have also shown that the models considered here 
are weakly coupled throughout the relevant parts of the 
parameter space, unlike some other models in the liter- 
ature that couple dark matter and dark energy and are 
ruled out by the quantum mechanical strong coupling. 

It is a testament to the many diverse sources of data 
in modern cosmology that the simple possibility of cou- 
plings between, say, dark matter and dark energy, can be 
constrained in so many different ways. The web of con- 
straints that we have delineated in this paper sets strict 



limits on allowed interactions in the dark sector, and may 
have important ramifications both for phenomenological 
models, and for fundamental theory. 
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